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ABSTRACT 

Context. The identification of pulsation modes in 6 Scuti stars is mandatory to constrain the theoretical stellar models. The non-adiabatic 

observables used in the photometric identification methods depend, however, on convection modeling in the external layers. 

Aims. To determine how the treatment of convection in the atmospheric and sub-atmospheric layers affects the mode identification, and what 

information about the thermal structure of the external layers can be obtained from amplitude ratios and phase lags in Stromgren photometric 

bands. 

Methods. We derive non-adiabatic parameters for 6 Scuti stars by using, for the first time, stellar models with the same treatment of 
convection in the interior and in the atmosphere. We compute classical non-gray mixing length models, and as well non-gray "Full Spectrum of 
Turbulence" models. Furthermore, we compute the photometric amplitudes and phases of pulsation by using the colors and the limb-darkening 
coefficents as derived from the same atmosphere models used in the stellar modeling. 

Results. We show that the non-adiabatic phase-lag is mainly sensitive to the thermal gradients in the external layers, (and hence to the 
treatment of convection), and that this sensitivity is also clearly reflected in the multi-color photometric phase differences. 
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•fh 1. Introduction 

%4 Asteroseismology uses the frequency of stellar pulsations as 
5^ ! a probe of the interior stellar structure, by comparing the ob- 
served frequencies with those determined from stellar models. 
In order to constrain the physical parameters of a star from its 
pulsation frequencies, it is necessary to identify the modes of 
pulsation, that is, to determine the spherical harmonic degree 
{, the azimuthal order m, and the radial order n for each mode. 
This is far from trivial in 5 Scuti stars, as their excited modes 
do not fall in the asymptotic range and, furthermore, rotation 
and possible coupling destroy any pattern regularity. As a con- 
sequence, mode identification is difficult and other complemen- 
tary observables must be used for this purpose. 

The mode identification photometric methods are based 
on the analysis of luminosity variations in different photo- 
metric passbands, and on the dependence of their amplitudes 
and phases on the spherical harmonic degree, I. These tech- 



niques are based on the theoretical study by Dziembowski 
d 1977b of the light variations generated by non-radial os- 
cillation. Afterwards, s e veral improvements were made by 
iBalona & Stobid dl979t). IStamford & Watsonl ill 98 lb IWatsonl 
d 1988b . iGarrido et al.1 dl990l) and iHevnderickx et all dl994h . 
The main shortcoming of the above mentioned approaches is 



that the non-adiabatic amplitude and phase of the effecti v e tem- 
perature are treated as free parameters. ICugier et al.1 ([1994) 



updated the method by using Dziembowski's non-adiabatic 
code to derive amplitude and phases and applied it to f3 Cep 
stars. This last approach did not take into account the non- 
adiabatic iegenfunctions in the atmosphere layers, and assumes 
the Lagrangian temperature variation to be equal to the local ef- 



fective temperature variation at the photosphere. Dup ret et al 



( 2002) showed, however, that the temperature perturbation in 
the photosphere is rather different from the variation of the 
local effective temperature, because of the significant optical 
depth variation produced by the oscillation. 

Since the non-adiabatic observables (multi-color magni- 
tude variation) strongly depend on the variation of the at- 
mospheric thermal structure, their accurate determination re- 
quires to solve the adequat e non-adiaba t ic equ ations in the 



stellar atmosphere as well. Dup ret et al.l (12003b included an 
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improved treatment of the oscillations in the non-grey atmo- 
spher e. The application of their model to f3 Cephei and SPB 
stars JDupret et al.ll2003l) an improvement in mode identifica- 
tion. 

Stars in the Hertzprung-Russell Diagram (HRD) region 
where 6 Scuti stars are located should have, according to theo- 
retical stellar models, an external structure with one or two con- 
vective regions due to H and He ionization. Since these regions 
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are very shallow, they are as well very over-adiabatic. As a con- 
sequence, the thermal structure of the outer layers, and there- 
fore the photometric amplitude ratios and phases differences in 
different bands, ar e very sensitive to the d etails of convection 
modeling, such as Balo na & Eversl dl999b showed in the first 
non-adiabatic study of 5 Scuti stars using the Dziembowski's 
code. 

The "standard model" of convection in stellar evolution is 



the mixing length theory (MLT Boh m- Vitensd 1 1 9581) . where 
turbulence is described by a relatively simple model. This con- 
tains essentially one adjustable parameter, a, relating the mix- 
ing length to the local pressure scale height. Several recent pa- 
pers have disc ussed the effect on the non-adiabatic observables 
of varying a dDaszvnska-Daszkiewicz et alj|2003t Mova et al 



2004 ; Daszvnska-Daszkiewicz et all 2005b . However, in those 
papers, the convection treatment is changed only in the regions 
computed with the corresponding stellar structure code, and not 
in the atmospheric layerfl where the multi-color magnitude 
variations due to oscillations are generated. 

The analysis of the effect of convection modeling on stellar 
structure and evolution is not always obvious, since the avail- 
able atmosphere models are tipically computed for a given the- 
ory of conv e ction and for a fixed value of c-MLT. Luckily, 



Heiter et al. 



(2002) have published new atmosphere model 
grids based on Kurucz's ATLAS9 code. In these new grid 
(thereafter NEMO grid), the atmosphere models have been 
computed for different treatments of co nvection: MLT (with 
two dif ferent values of the a parameter). ICanuto & Mazzitellil 
dl991al) (CM), and ICanuto et al.1 d 19961) (CGM) treatments 
(also known as FST - Full Spectrum of Turbulence- treat- 
ments). The s e choi ces are extensively motivated in Sect. 2 of 



Heiter et al. 



(2002), and the main characteristics and differ- 
ences among these convection treatments will be described in 
Sect. 

By including the NEMO grids in our model computations, 
we are able to study in a self-consistent way the effect of con- 
vection treatment on the oscillation amplitudes and phases in 
different passbands. The analysis presented here is consistent 
at two levels: first, the stellar models (interior and atmosphere) 
used in the non-adiabatic computations have been obtained by 
using the same (either FST or MLT) convection treatment; and 
second, the color transformations tables, as well as the limb- 
darkening law derivation, are based on the same atmosphere 
models as those we included in the stellar structure modeling. 
In Sect. 12 we present the models used in this study: the inte- 
rior equilibrium models (Sect. l2~TT i. the atmosphere models and 
limb-darkening (Sect. 12.3b . the non-adiabatic oscillation code 
(Sect. 12.4b and the color transformations giving the observable 
properties of S Scuti pulsations (Sect. 12.5b . In Sect. 12.21 we 
provide an outline of the problem of the modeling of convec- 
tion in A-type stars. 



One of the main results of iHeiter et al.l (12002b is that MLT 
models with a small mixing length parameter (e.g., a ~ 0.5) 
are equivalent, in the atmospheric region where the observed 
flux originates, to FST models. Both treatments predict, in- 



deed, a low convective efficiency for these layers. The ther- 
mal structure in the subphotospheric layers for models with 
MLT{a = 0.5) and FST may be, however, quite different. We 
will study in this paper if the non-adiabatic features of 5-Scuti 
oscillations are sensitive to these differences. 

In Sect. [3] we compare the results obtained with MLT and 
FST treatments and explain the origin of the differences. We 
analyze as well the respective weights of the over-adiabatic 
structure and of the color transformations on the theoretical 
amplitude ratios and on phase differences in Stromgren pho- 
tometric passbands. Finally, we consider in Sect. H] the appli- 
cation to a particular 5 Scuti star: AN Lyn, and, in Sect. [5] we 
present our conclusions. 



2. Theoretical models 

2.1. Stellar models 

The stellar models were computed by means of the stellar evo- 
lution code ATON3 .0, recently updated to be used in asteroseis- 
mology modeling ( D'Antona et al.ll2005l) . We computed mod- 
els for a metal mass fraction Z=0.02 and assumed a helium 
mass fraction Y = 0.28. At temper atures T > 12000 K we a dopt 
the OPAL radiative opacities (k) (Iglesias & Rogers 1996|) for 
the solar Z-distribution from lGrevesse & Noelsl dl993b . In the 
high-density (p) regions the opacities are linearly extrapo- 
lated (log k vs. l og p), and harmonicall y added to conduc- 
tive opacities by Itoh & Kohvama (1993). At lower tempera- 
tures we use lAlexander & Ferguson! d 19941) molecular opaci- 
ties (plus electron conduction in full ionization) for the same 
H/He ratios as in the OPAL case. A complete description of 
th e equation of state (E OS) used in ATON3.0 code is given 
in 



D'Antona et al. (2005). In short th is code uses OPAL( 2001) 



EoS (Rog ers & Navfonovl2002"l) and lSaumon et alJd 1995b EoS 
for low temperature and high density domain. 

Several treatments of convection are available in ATON3.0: 
the classic m ixing length theory (MLT) model, with the 
Cox & Giulil (1968) for mulation, and two d ifferent 
models, according to Canuto & Mazzitelli ( 1991bb 



FST 
and 

Canuto et al.l (11996) approaches. A detailed des cription of the 



convec tion modeling in ATON3.0 is provided by Ventura et al 



; 1998). For the present study, we adopt MLT and ICanuto et al 



19961 thereafter FST or CGM) models. Except for explicitly 
mentioned cases, the a-MLT and the a* -FST parameters are 
the same as those in the atmosphere models which were used 
to derive external boundary conditions. 

As surface bou ndary conditions (for detail, see 



Montalba n et all |2004) we take the NEMO atmosphere 



1 The exact optical depth, r = 2/3 or t = 10 depends on the partic- 
ular stellar evolution code used 



grids. As mentioned above, these new atmosphere models 
have been computed using MLT model (with a — 0.5 and, for 
a smaller range of parameters, a - 1 .25 as in cla ssic ATLAS9 
mode l), and, as well, CGM with a* = 0.09 dCanuto et al. 
Il996b . That allowed us to consistently build stellar models 
(interiors and atmospheres). The boundary is located at optical 
depth T p h = 10, as suggested by Hei ter et al.l (2002)), except 
when MLT models are computed with a or-MLT value in the 
interior different from that of the atmosphere. In this case, the 
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boundary is located at T p h = 1, to avoid mixing of different 
convection efficiencies in the overadiabatic region. 

Helium and metal gravitational settling are not included in 
these computations. All the models are computed following the 
evolution from the pre-main sequence (PMS) phase. As nu- 
merical resolution is concerned, the number of models for the 
PMS phase is of the order of 600, and the evolution from the 
ZAMS to the TAMS is done in -350 models, besides, 2500- 
3000 mesh points are used to describe the stellar structure. 

Before applying non-adiabatic computations, the stellar 
structure of the equilibrium model is extended by matching the 
atmosphere structure (from r p h up to t — 10~ 4 ) as obtained by 
interpolation in the grid of atmosphere models. 

2.2. Convection modeling in A-type stars 

Standard models for A-type stars predict stellar structures with 
a convective core, and one or two external convective regions, 
the inner one corresponding to Hell ionization, and the shal- 
lower and outer one to that of HI and Hel. Convection in the 
central region is almost adiabatic, while the external convective 
regions are highly overadiabatic, but very thin and shallow. As 
a consequence, the details of convection in A-type stars (except 
for the overshooting that changes the dimension of the convec- 
tive core) have a negligible effect on the their location in the 
HR diagram. 

Nevertheless, the convective heat transport changes the 
temperature gradient, and affects, therefore, the spectral line 
profiles and the emission in different photometric bands. In 
fact, there are spectroscopic evidences of the presence of con- 
vection in A-type stars, such as the peculiar behavior of Balmer 
lines as a function of temperatu re, the sh ape of bisectors, or 
chromospheric activity (see lSmalleyl2004l for a review). These 
observations suggest that convective velocity fields of several 
km/s are present in A-stars, while the temperature gradient 
must be close to the radiative one, except for the coolest ob- 
jects (~ 7000 K). 



It has been already s hown by [ Fuhrmann et al. (1993), 
Fuhrmann eU d. ( 19941). van't Veer-Menneret & Megessier l 



dl996|) . Ivan't Veer-Menneret et al J dl998h. iBarklem et al.1 
hood). ISmallev & Kupkal dl997t) . Ismallev & Kupkal ( 1 1998b . 



Smallev et al] (12002b . that spectroscopic and photometric fea- 



tures seem better fitted by adopting low efficiency convective 
models (for instance, MLT with o , =0.5 or CGM) which provide 
temperature gradient close to the radiative one. Nevertheless, 
neither MLT or CGM succeed in matching all the photometric 
features in the A-F stars effective temperature domain. 

Other observational evidences in favor of low 
efficient con vection in A-type stars h as been pro- 
vided by 



Daszvhska-Daszkiewicz et al 
Daszvhska-Daszkiewicz et alj ( 2005b in 



(120031) and 
their fits of the 



seismic features of several 6 Scuti stars. 

Either, non-loca l models of convection 

( Kupka & Montgomery 2002[) or 2D and 3D numerica l 
simulations dLudwig et alj 1 19991: iFrevtag & Steffenl |2004|) 
indicate that, in order to reproduce the convective flux in the 
external layers, the value of the MLT-a parameter should 



vary from 0.5 at 8000 K to ~1 or larger at 7000 K. On the 
other hand, the a value required to reproduce the flux in the 
Hell convective zone should also be four times larger than 
that used in the HI, Hel convective zone. For tunately (or 



unfort unately) for what concerns the present study. [Mova et al 



(2004) has shown that a change of convection efficiency in the 
Hell convective region does not affect the excited modes of 
6 Scuti stars. 

In the computation presented in this paper we use the 
Mixing Length Theory (Bohm-Vitense 1958), and a the "Full 
Spectrum of Turbulence" formalism (CGM). The main im- 
provement of CGM with respect to MLT is that while the latter 
treats the heat transport by mimicking the spectral distribution 
of eddies by one "average" eddy (reliable only for high viscos- 
ity fluids), the former overcomes the one-eddy approximation 
by using a turbulence mode l to compute the full spectrum of 
a turbulent convective flow. Canuto (1996) showed that a con- 
sequence of this different heat transport is that in the limit of 
highly efficient convection (S » 1, where S is the convec- 
tive efficiency, see e.g. CM for details) MLT underestimates 
the convective flux, while - in the low efficiency limit - MLT 
overestimates it. So, CGM convective fluxes are ~ 10 times 
larger than the MLT ones for the S » 1 limit while, in the low 
efficiency limit (S «C 1), the CGM flux is ~ 0.3 times the MLT 
one. This behavior yields, in the over-adiabatic region at the 
top of a convection zone, steeper temperature gradients for FST 
than for solar-tuned MLT models. In inefficient convection, the 
temperature gradient sticks to the radiative one and starts split- 
ting only when convection becomes efficient. Another differ- 
ence between CGM and MLT concerns the length scale (A) for 
transport processes: Amlt - a ■ H v and Acgm = Z + a*H p i op , 
where z is the distance to a boundary of convective region, and 
the second term is a fine tuning parameter that allows small 
adjustments, if e xact stellar radii are needed ( e.g., in helio- 
seismology). As ICanuto et al. I dl996b stress, the role of a* in 



solar FST models is radically different from that of a in the 
MLT model. In fact, FST tuning affects only layers close to 
the boundaries. For inner layers, z quickly grows and becomes 
much larger than a*H Pttop . In the case of A-type stars, however, 
the small dimension of convective zones implies that this term 
could be a significant fraction of z. For the models presented in 



this paper, we use a* = 0.09, as derived by Canuto et al. ( 1996) 



from solar calibration^ and adopted by Heiter et al. (2002) for 
their atmosphere models. 

2.3. Atmosphere models and limb-darkening 

The original ATLAS9 code treats convection by using MLT 
with a = 1 .25 a nd thas the possibili ty of including a sort of 
overshooting (see lCastelli et aljl997alb[ for details) The value 
a=1.25 and the option of "approximate overshooting" were 
originally adopted to fit the intensity spectrum a t the center of 



the solar disk and the solar irradiance. However, Castelli et al 



dl997bb showed that these quantities are much more sensi- 
tive to the overshooting-on mode than to the value of a it- 



2 Note that this value was derived for the solar calibration using 
grey-atmosphere as external boundary conditions 
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self (see also iHeiter et alj|2002l). T he new grids of ATLAS9 
atmospheres by Heiter et al.l ( 2002 ), which we will use in our 
computations, introduce som e improvemen ts w ith respect to 
the or iginal ones published bv lKuruczl (1 1993b and Castelli et al 



( 1997b), in particular, they provide a finer grid spacing (AT e s, 
Alogg), as well as a higher vertical resolution. For the reader's 
convenience, we recall the main characteristics of these grids: 



i) Three different models of convection are considered: 
MLT(a = 0.5), FST according to CM and to CGM. 

ii) For FST models, the vertical resolution in the atmospheric 
integration is increased from 72 to 288 layers, spanning the 
range from log tr oss = —6.875 to log TR OSS = 2.094 (where 
tross is the average Rosseland optical depth). 

iii) Effective temperature range: from 4000 to 10000 K, with 
AT eS = 200 K, and 

iv) gravity ( logg) from 2.0 to 5.0, with Alogg = 0.2. In 
the n ew ATLAS9-ML T gri ds, the convec t ion is described 
as in 
V/A = 



Castellil d 19961) and ICastelli et al] dl997bh (that is 
A/6 and y = 0.5), but with a = 0.5 and without 
"overshooting". 

In the C GM model atmospher es, the convective flux is com- 
puted as in ICanuto et al. (1996), but the characteristic scale 
length is defined as: A = min(z top + a* H Pttop , z bo t + a*# P ,bot) 
where the index "top" and "bot" refers to top and bottom of the 
convective region, and a* = 0.09 (see footnote 2). 

The photometric methods (Sec. 12.51 ) used for 6 Scuti mode 
identification are based on an analytical expression of the 
monochromatic oscillation amplitudes and phases. In this ex- 
pression (Eq.Q]) explicitly appear the weighted limb-darkening 
integrals {bcj, Eq. [2]) and their derivatives with respect to 
log r e ff and \ogg (Eqs. [3]and|4]i. The new and denser grid of 
model atmospheres allows to derived smoother bcj, and their 
derivatives as required for mode identification method (Garrido 
2000). 

In this paper we will use the limb-darkening coefficients 
(LDC) for th e Stromgren photometric system, as derived by 



Barb an et al.l (2003) by using the quadratic-law, for a sam- 
ple of new atmosphere models with effective temperature be- 
tween 6000 and 8500 K, \ogg between 2.5 and 4.5, and so- 
lar metallicity. These authors determined, for each uvby band, 
the best limb-darkening law as well as the integrals bcj and 
found similar results for CGM and MET (ff=0.5), but signifi- 
cant differences between CGM and the usual MET (or= 1 .25). 
Besides, they showed that the differences in the LDC's fol- 
low those in the vertical temperature structure of the atmo- 
sphere model due to the different treatments of convective 
transport. The differences between CGM and MET (cr=1.25) 
models are higher for low temperatures, but the differences 
propagate toward high temperatures for model atmosphere with 
higher gravities. The differences between LDC's computed for 
MLT(ar=0.5) and CGM atmosphere models are much smaller 
than for the MLT(a:=1.25), but they can reach the same order 
of magnitude as for MLT(a-l.25) when intermediate gravity 
and low temperature (~ 7600 K) are considered. As a gen- 
eral statement, the differences in the bf coefficients and deriva- 
tives computed with CGM and MET (ff=0.5) are much smaller 



(~ 5%) than differences obtained between computations with 
CGM and MLT (a=1.25) (~ 20%). 

As for the LDC's the effect of convection treatment on be 
is maximum at low effective temperature and low gravity mod- 
els, and these large differences extend towards higher tempera- 
tures for models with high gravity. Furthermore, the magnitude 
of the differences in be due to the convection treatment is I— 
dependent. 

The use of a higher resolution in optical depth in the model 
atmospheres and the different treatment of the convection im- 
ply significant changes in the LDC values. In this paper we use 
these results to analyze the effect of LDC's and of the external 
thermal stratification on the theoretically determined amplitude 
ratios and phases. 



2.4. Non-adiabatic treatment 

We do a non-adi abatic seismic anal ysis of our models by using 
the code MAD dDupret et al.ll2003l) . In this code the adiabatic 
iegenfunctions are used as trial input for the non-adiabatic com- 
putations. The input model is a stellar model extended with the 
atmosphere corresponding to its 7^, logg and adopted con- 
vection treatment. In this context, we use for the first time the 
Heiter et al.l (120021) atmosphere models. 



Since the convection-pulsation interaction implemented in 
MAD cannot be applied to the FST formalism, we adopt in our 
study the frozen convection flux approximation, that is, we ne- 
glect the Lagrangian variation of the convective luminosity (ra- 
dial component) and the Lagrangian variation of the transver- 
sal component of the convective flux. This approximation is not 
appropriate to describe the red edge of 6 Scuti instability strip. 
Nevertheless, given the comparative character of this work and 
the low efficiency of convection of most of the models consid- 
ered here 0, we expect the convection-pulsation interaction to 
be not relevant to this study. The transfer equation treatment in 
the MAD code is differen t in th e interior and in the pulsating 
atmosphere (Dupr et et al. 20031) . In the latter, the hypothesis 
at the base of the procedure is not justified for convective re- 
gions, and therefore, the match layer between the interior and 
the atmosphere must be located outside the convective enve- 
lope. Even if in the extended stellar model the match between 
atmosphere and interior is done in a convective region, and the 
properties of a given layer come formally from the atmosphere 
model, the non-adiabatic treatment reserved to the atmosphere 
will be apply only to layers such that the convective to radiative 
luminosity ratio is smaller than 10~ 10 . 

These computations make possible to derive the phase lag 
(tffj) between the local relative variation of the effective tem- 
perature and the relative radial displacement (6 r) and the cor- 
responding amplitude ratio (/x), which are basic ingredients 
to determine the magnitude variations in different photometric 
passbands (see next section). 



FST models and MLT ones with a parameter smaller than one. 
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2.5. Multi-color photometry 
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Fig. 1. Fraction of luminosity transported by radiation with re- 
spect to total luminosity as a function of temperature for sev- 
eral 1.8 M Q models in the same HRD location (7^=7400 K, 
logg=3.98) but computed with different treatments of convec- 
tion. The FST model has been compu ted with FST conve ction 
in the interior and FST atmospheres by Heiter et al.1 d2002l) . The 
MLT(a; = 0.5) mo del implies a-0-5 in the interior and the cor- 
responding g=0.5 iHeiter et all (2002) atmosphere. The MLT 
models with a + 0.5 were computed with the indicated a value 
in the interior up to the optical depth r = 1 and MLT(o' = 0.5) 
atmosphere model for r < 1 . 



The non-adiabatic quantities iffj - (piST^/T^) - 4>(6r/R), 
fr = \8TeslTtfs\l\6rlR\ and f % = \Sg eS /g eS \ (where g eff is the 
effective gravity and 5g the variation of gravity at the pho- 
tosphere) can be related to photometric observables such as 
the amplitude ratios and phase differences between magni- 
tude variations in di fferent photom etric passbands. It has been 
shown, in fact, (e.g. Watson 1988, and references therein) that, 
in a one-layer linear approximation, the magnitude variation at 
a wavelength A produced by a stellar pulsation with spherical 
harmonic degree t, azimuthal order m and angular frequency <r 
is given by 



5mx 



2.5 



a P';'(cos b a { - (( - 1) (€ + 2) cos (erf) 



In 10 

+ / T cos (at + fx) (a TA + Pta) 
- f g cos (at) (a^ +/3 gA )} 



(1) 



where i is the inclination angle between the stellar rotation axis 
and the observer line of sight; a is the amplitude of the relative 
radial displacement at the photosphere and P, is the associated 
Legendre function of degree i and azimuthal order m. ba is 
related to the limb-darkening law by 



f 

Jo 



(2) 



where hx is the normalized limb-darkening function, and fi 
cos with 9 the angle between the line of sight and the normal 
to the local stellar surface. The rest of quantities appearing in 
Eq. [Tj ajA, /3ta, a g A and /3 g A are derived from the appropriate 
atmosphere model with effective temperature r e jf gravity log g 
and given convection treatment: 
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Fig. 2. Temperature gradient as a function of temperature inside 
the star, for some of the models in Fig.Q] 
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An appropriate way to test multi-color theoretical predic- 
tions is to construct phase-amplitude diagrams correspond- 
ing to well-chosen combinations of photometric bands (e.g. 
Garrid o et al. 1990). The theoretical results corresponding to 
modes of different spherical degrees t occupy well distinct re- 
gions. This makes possible the identification of I by search- 
ing for the best fit between theory and observations. The non- 
adiabatic quantities f T and i^ T , which play a mayor role in 
Eq.Q] are highly sensitive to the characteristic of t he convective 



envelo pe (see. Fig. [7] and [8). On the other hand. iBarban et al 
(2003) showed that the sensitivity of ba and its derivatives 



to the convection model is different for different photometric 
passbands and different degree I. We expect, therefore, that 
constraints on the characteristics of the convective envelope 
could be obtained from very precise multi-color photometric 
observations. 
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3. General results for 6 Scuti stars 

3.1. Equilibrium models 

As already mentioned, thanks to the shallow convection enve- 
lope in A-type, the HRD location of a stellar model of fixed 
mass and chemical composition is not affected by the treatment 
of non-adiabatic convection layers. This fact allows us to easily 
separate the effect of convection. We have computed models for 
1.8 and 2.0 M with FST formalism, and with MLT(a = 0.5). 
We have also computed MLT models with a ml =0.1, 0.3, 0.7 
and 1.7fl Even if for a given central hydrogen content, all 
the computations yield the same point in the HRD, the char- 
acteristic of the external layers are affected by different con- 
vection treatments. The two main physical quantities affected 
by this treatment are the energy fraction transported by con- 
vection (Fig. [TJ and the temperature gradient (Fig. |2}. As it 
is well known, the energy fraction transported by convection 
increases as a increases, while the temperatu re grad i ent de - 
creases. On the other hand, as pointed out by ICanutol dl99q) . 
the energy fraction transported by convection in FST models 
is much smaller than that of MLT treatment in the case of low 
efficiency convection. 

In next section we analyze the consequences of these dif- 
ferent structures on the non-adiabatic observables. 

3.2. Non-adiabatic results 

We present now the non-adiabatic results obtained for mod- 
els with a given mass ( 1.8 M Q ) and effective temperature 
~7400 K. In Figs. [3] and [4] we plot the relative amplitudes 
(/t) and phases (tf/ T ) as a function of the pulsation constant 
Q.Q = P ^(Rq/R)\M/Mq), with P being the period in days. 
In 5 Scuti stars Q 0.033 days for the fundamental radial 
mode. The models used in Figs. [3] and |4] are the same as in 
Figs. Q] and [2] for the same kind of convection, for all of them 
Tdyn/r dyn ,o = ^(R/Rq) 3 (Mq/M) = 3.104282. 

The results presented in Figs.[3]and|4]confirm that the non- 
adiabatic quantities (particularly the phase-lag xf/i) are very 
sensitive to the convection treatment. As the value of the MLT- 
a parameter increases, the phase-lags decreases. Furthermore, 
for low values of the MLT-a' parameter, the i^t value for high 
radial order modes rem ain close to t he adi abati c one (180°), 



as alre ady discussed by iMova et alj (120041) and iDupret et al 
d2005l) . 



fj and iffi values for FST model are in between those ob- 
tained for MLT(a = 0.1) and MLT(a = 0.3) models. In all our 
models except one, we use NEMO atmosphere models with 
a low efficient convection (FST or MLT-a=0.5)). The excep- 
tion is for the model with an interior structure computed with 



4 We would like to call attention on the fact that since the available 
atmosphere models use MLT(o'=0.5) or CGM, we are able to com- 
pute only "complete" MLT(or=0.5) and "complete" FST models (the 
same treatment of convection in the interior structure and in the atmo- 
sphere). The internal structure of models labeled with MLT a=0. 1, 
0.3, 0.7 and 1.75 has been computed until r = 1 with these a 
values, and with boundary conditions for temperature and pressure 
((T, P){t =1)) provided by the the MLT(a;=0.5) atmosphere model. 
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Fig. 3. The non-adiabatic quantity fj = \6T e ff/T e g\ as a func- 
tion of the pulsation constant Q (in days) for modes with spher- 
ical degree I = 0, 1, 2 and 3. These modes correspond to 1 .8 M 
models in the same HRD location and computed with different 
convection treatment. The meaning of labels is the same as in 
Fig. [□ The label MLT (ar=0.5) Kurucz (1.25) refers to a model 
computed wit h MLT(ar=0.5) in the interior computation and 
Kuruczj(11993|)'s atmosphere models. 



MLT(a = 0.5) up to t = 1, and with the boundary conditions 
and the atmosphere structure given by a — 1 .25 Kurucz's atmo- 
sphere models. The results provided by this last model and by 
the "complete" MLT(ar = 0.5) model are significantly different. 
This shows the relevance of treating consistently the interior 
and the atmosphere. 

We remark that the pulsation frequencies are almost iden- 
tical for the different convection treatments. The reason is that 
these models only differ in the outermost layers which have no 
weight at all on the frequencies of low order p-g modes. As 
the effect of temperature gradient in mode excitation is con- 
cerned, we state that for models with very inefficient convec- 
tion, i.e. FST, as well as MLT(or=0.1) and MLT(or=0.3), only 
modes with Q > 0.014 d are excited, whereas for MLT(ff=0.5) 
this limit moves to higher frequencies, and modes with Q > 
0.0125 d are excited as well. 

In order to understand the behavior of iffj for different 
convection models, we plot in Fig. [5] the phase-lag between 
the total luminosity perturbation and the radial displacement 
(<p{6L) - <p{5r)) for the radial fundamental mode in the same 
models as in Figs. [3] and |4] In deep regions (log T > 5) the 
pulsation is quasi-adiabatic and the phase-lag is 180°. A first 
phase-lag, quasi-independent of the convection treatment, oc- 
curs at log T = 4.8^.6, in the He n partial ionization zone, 
and a second one appears in the H partial ionization zone 
(4.2 > log T > 4.). This second phase-lag increases as the MLT 
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Fig. 4. As in Fig. [3] for the phase-lag (in degrees) between the 
local effective temperature variation and the radial displace- 
ment. 

a parameter increases, which leads to phase-lag values close to 
0° at the stellar surface. 



As shown in Dupr et et al.l d2005l) . the phase-lags can be in- 
terpreted on the basis of the following integral expression: 



J [8L/L j L 

J [6L/L J LdlnT/dr 



(5) 



This equation makes evident the different contributions to 
the phase-lag through the star. Thus, a significant contribution 
can be generated in a region where the thermal-relaxation time 
is of the same order or larger than the pulsation period, and/or 
in a region where the pulsation is highly non-adiabatic. This 
last situation is found in the convective zone coinciding with 
the H partial ionization zone (HCZ), where %{6s L/{c v 6L)} has 
a high value due to the the large opacity bump and to the tem- 
perature gradient. The sensitivity of the phase-lag to the con- 
vection treatment can be made more obvious by noting that the 
iegenfunctions 6L/L when written as a function of T instead of 
as a function of mass, are essentially the same for different con- 
vection treatments in the equilibrium models. The main reason 
is that the opacity is essentially a function of temperature and 
is not sensitive to a for a given interv al or temperature (see also 
Mova et al.ll2003l:lDupret et al.l2005b . As shown in Fig.|H what 
significantly changes from one model to another is the temper- 
ature gradient. When a increases, the temperature gradient in 
the HCZ decreases quickly. According to Eq. [5] (second line) 
the phase-lag contribution from this zone increases, therefore, 
with a, and hence, the phase-lag can be considered as an indi- 
cator of the temperature gradient in the HCZ. 
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Fig. 5. Luminosity phase-lag (4>{5L/ L) - <p(6r/r) (in degrees) as 
a function of temperature inside the star for the radial p\ mode 
corresponding to the same 1.8 M models than in previous fig- 
ures. 
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Fig. 6. As Fig.|5]for the n=4 overtone. 



In Fig. 12 we see that the mean temperature gradient of the 
FST model is between the temperature gradients of the MLT 
models with a = 0.1 and a = 0.3, and consequently the phase- 
lags predicted by FST models are between the a = 0.1 and 
a = 0.3 MLT ones (see Figs|5]and |6]). 

It is well known that different modes probe the different 
layers of stellar interiors. Similarly we can expect, in principle, 
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that the non-adiabatic phase-lags predicted for different modes 
probe the temperature gradient at different depths. It is evident 
from Fig.|2]that the exact shape of the temperature gradient pre- 
dicted by FST models cannot be reproduced by MLT models, 
whatever the value of a. Consequently, the phase-lag predicted 
by FST models cannot be reproduced by the same MLT model 
in the whole frequency spectrum. This can be seen in Fig. [4] For 
periods close to that of the fundamental radial mode, the FST 
phase-lags are close to the MLT(a = 0.3) ones. But for smaller 
periods, the FST phase-lags are closer to the a = 0.1 MLT ones. 
In Fig. [6]we plot <p{6L) - <p(6r) of the p4 modes computed for 
the same models of Fig. [5] As in the case of the pi mode, we 
have d<p(5L)/dr < in the partial ionization zones of He and H, 
but, at variance with it, d<p(6L)/dr > in the intermediate re- 
gion between the He and H partial ionization zones, so that the 
phase-lag of the p4 mode is closer to 180°. These positive val- 
ues of d<p(6L)/dr are associated with the damping of the modes 
in this region. Different shapes of the iegenfunctions, depend- 
ing on the mode, imply a different sensitivity of the phase-lag 
to the convection treatment. Nevertheless, the differences be- 
tween the phase-lags of the FST and the a = 0.1, a — 0.3 MLT 
models are small whatever the mode, and it would be too opti- 
mistic to conclude that the present precision of observations for 
the phase-lags would make possible to discriminate between 
FST and MLT models with the lowest a values (see discussion 
below). 

Finally, in FigslTHTOl we present the sensitivity of fj and 
i//t to the convection treatment as a function of the effec- 
tive temperature, for 1.8 and 2.0 M main-sequence models. 
Comparing the MLT and FST curves, we see that the FST re- 
sults are always between the a = 0.3 and a — 0.1 MLT re- 
sults, closer to a = 0.1 for hotter models and to a = 0.3 
for cold models. It is also evident that is much more sen- 
sitive than /x to the convection treatment in the equilibrium 
model. Furthermore, the phase-lags decrease as effective tem- 
perature decrease, and tend to zero more rapidly as convec- 
tion efficien cy increases. The latter re sult agrees with th ose 
obtained by iMoya et"al] d2004 and by bupret et all d2005l) in 
the frozen convection approximation. We remark that if/j differ- 
ences among different convection treatments are larger for the 
cold models than for the hottest ones. These results can be un- 
derstood as follows: when the effective temperature decreases, 
the size of the convective envelope increases quickly and the 
gradient of temperature decreases. This is similar to increasing 
a and, for the same reasons, the phase change between bottom 
and top of convection zone increases. A second effect of de- 
creasing effective temperature is to increase the opacity bump. 
Consequently, the associated bump of the non-adiabatic iegen- 
functions is larger and implies as well larger phase changes. 

if/j for FST models follows the MLT(q'=0.1) behavior un- 
til r e ff~ 6800 K, then it changes suddenly its slope to follow 
the MLT(ar=0.3) <At curve. This is explained by the properties 
of FST convection (see e.g. iD'Antona et alj|2002l) which in- 
duce a rapid change of the depth of the convective zone at that 
J e ff. That change from shallow to deep convection has the same 
effect on tf/j as the increasing of a parameter in the MLT mod- 
els. We note that FST if/j values for r e ff > 6500 K models are 
between 90 and 140°, in good agreement with the typical ob- 




6500 7000 7500 8000 6500 7000 7500 8000 



Fig. 7. The non-adiabatic quantity fj = \6T e ff/T e fi\ for four dif- 
ferent modes (■£ = 0, n = 1; ( = l,n = \;l = 2, n = and 
n = 1) as a function of the effective temperature for the MS 
evolution of 1.8 M models computed with different convec- 
tion treatments. The meaning of the labels is the same as in 
Fig.Q] 



servational values in 6 Scuti stars, and with the th eoretical pre 



diction s of time-dependent-convection models by Dupret et al 
d2005l) . 



For 2.0 M Q models, where the convective zone is shal- 
lower than in 1.8 M Q ones, the effect of considering FST or 
MLT(or=0.5) is smaller than for 1.8 M models, for both the 
amplitude ratio and the phase-lag. Nevertheless, FST models 
show always iffj > 100° while for MLT(o:=0.5) decreases 
down to ~ 50° for lower effective temperatures. High a values 
provide as before, too low phase-lags. 

3.3. Multi-color photometry 

The above mentioned results directly affect the phase differ- 
ences and amplitude ratios between different photometric pass- 
bands. More precisely, as introduced in Sect. 12.51 two differ- 
ent ingredients allow us to determine the magnitude variations 
in different passbands: the non-adiabatic normalized ampli- 
tudes and phases of effective temperature variation (fj and i/r-r); 
and the dependence of the monochromatic flux and of limb- 
darkening on the effective temperature and gravity. All these 
quantities, as mentioned in the previous sections, depend on 
the adopted treatment of convection. 

We con si der h ere three families of atmosphere mod- 
els: iKuruczl dl993l) 's models with MLT(g-1.2 5 ) trea tment, 
and NEMO atmosphere models dHeiter et al.l 120021) with 
MLT(a=0.5) and FST. 
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Fig. 8. As Fig|7]for the temperature phase-lag: iffj. The gray 
bands indicate the models for which the corresponding l,n 
mode is excited. 



Fig. 10. As Fig.[9]for the temperature phase-lag: ijjj. The gray 
bands indicate the models for which the corresponding l,n 
mode is excited. 
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Fig. 9. As Fig.|7]for 2.0 M models and for only three different 
convection treatments. 



We first recall that, because of the dependence of Eq. ([T) 
on the degree I of the mode, the comparison between theo- 
retical and observed photometric amplitude ratios and phase 
differences makes ^-identification possible. On the other hand, 
different convection treatments lead as well to different results, 
so that it is important to check how this affects mode identifi- 
cation. The classical diagrams for mode identification plot the 



amplitude ratio, A(t,- y )/A y , versus 0(b- y > - 4>y Fig. [TT] shows 
the theoretical predictions for models with r e fj = 7400 K, for 
modes with degrees i — to 3 and frequencies close to the 
fundamental radial mode. The results correspond to different 
models: one with MLT(ar = 0.5) treatment in the interior and 
in the atmosphere, a second one with MLT(o' = 0.5) in the in- 
terior but MLT (a = 1.25) in the atmosphere and a third one 
with FST treatment both in the interior and the atmosphere. We 
see that the regions corresponding to different I are still clearly 
disjoint in this diagram when different convection treatments 
are considered. Moreover, we note that the typical precision 
of photometric methods, i.e. a few degrees for the phases and 
a few percents for the amplitude ratios, would allow in prin- 
ciple to discriminate between different convection treatments. 
For the highest frequencies, however, this discrimination will 
be possible only for I > 2 (Fig.[T2l. 

A second important point is to examine the respective ef- 
fects of the non-adiabatic calculations and of the atmosphere 
models on photometric amplitudes and phases. In Figs. [T3l and 
[T4l we plot, for three model families, the amplitude ratios and 
the phase differences between two passbands as a function 
of the effective temperature. Full circles correspond to mod- 
els where: the FST treatment of convection is adopted both in 
the interior and the atmosphere, and the monochromatic fluxes 
and limb-darkening, used in the transformation from theoreti- 
cal to observational plane, were computed with the same kind 
of model atmosphere as in the stellar modeling. Open squares 
correspond to "complete" MLT(o' = 0.5) models, and as in the 
previous case, to a consistent transformation to the observa- 
tional plane. Finally, the crosses correspond to a hybrid case 
where we take the non-adiabatic results fj and xf/j of "com- 
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Fig. 11. Color to light amplitude ratio vs the phase difference 
diagram for three 1.8 M models: complete FST model (cir- 
cles), compl ete MLT(g-0.5 ) model(squares), and MLT(o'=0.5) 
model with iKuruca dl993l) 's atmospheres (crosses). Different 
groups correspond to different (: I = 3, upper-left; ( — 2, 
lower-left; I = 1, lower- middle and t — 0, lower-right. 



plete" MLT(q' = 0.5) models, but we use the color transforma- 
tions (Eq. CO)) obtained with the FST atmospheres. 



Fig. 13. Photometric amplitude ratio (upper panels) and phase- 
lag (lower panels) for the u and v passbands for two differ- 
ent modes (€ = 0, n = 1 and I = 2, n = 0) vs Teg 
for the MS evolution of 1.8 M models. Solid circles corre- 
spond to complete FST models with the color transformation 
provided by FST atmosphere models. Open-squares refers to 
complete MLT(a=0.5) with the color transformation given by 
MLT(ff=0.5) atmospheres; and crosses corresponds to com- 
plete MLT(ar=0.5) and the color transformation given by FST 
atmosphere models. 
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Fig. 12. As Fig. Qj] but for a mode with pulsation constant 
Q=0.015 



Fig.[T3lshows the amplitude ratio and phase differences be- 
tween m and v passbands, and Fig. [14] those between b and y 
passbands. Note that scales and ranges in those figures are not 
the same and that the effect of different convection treatments 
are much larger for the u and v passbands than for the b and 
y ones. We remark in both sets of figures an other interest- 
ing result: MLT and hybrid models show very close phase dif- 
ferences, and hence, the photometric phase differences mainly 
provide information about the non-adiabatic results. As shown 
in the previous section, if/j essentially depends on the tempera- 
ture gradient in the HCZ. Therefore, this gradient can be con- 
strained by comparing the theoretical and observed multi-color 
photometric phase differences. 

The differences between MLT and FST treatments are 
greatest at r e ff~6800 K. We choose a 1.8 M Q star, with 
log r e ft=3.83 and logg=3.81, and we compare the effect of 
different external layers treatment on the monochromatic phase 
differences and amplitude ratios. In fig. [15] we show the phase 
differences for the ubvy bands as a function of the wavelength 
(a separate panel for each degree £ = 0, 1, 2 and 3). Different 
curves corresponds to "complete" FST and MCIXa^O.S) mod- 
els, MLT with FST color transformations, and MLT(q'=0.5) 
with Kurucz-MLT(a=1.25) atmospheres and the correspond- 
ing color transformations. 
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Fig. 14. As Fig. [13] but for the b and y passbands. 

We confirm that the atmosphere models used in color trans- 
formation do not have a significant impact on the monochro- 
matic phase difference (even if Kurucz-MLT(ar=1.25) tables - 
not included in Fig. [TBI - are used), though the effect increases 
with I. The same type of plots of fig. Q3] are displayed in 
fig. [16] for the monochromatic amplitude ratio. As it was al- 
ready shown in Figs. [3]and[4] the differences in fj and be- 
tween the different treatments of convection, decrease as mode 
frequency decrease. That is the reason why the ratios A x /A y and 
phase differences <p x - <p y for the modes with frequencies clos- 
est to the fundamental one presented in figs. [TBI and [TBI are so 
similar for the different treatments of the external layers; and, 
why the differences increase with the frequency (or degree £). 
We think that we must not over-interprete the results obtained 
for models mixing MLT(a=0.5) in the interior with Kurucz- 
MLT(o'=1.25) atmospheres. First, because, as already mention 
in Sect. 13.21 by mixing very different convection treatments, we 
get unforeseeable results such as those shown in Fig. |4] where 
the large value of is probably due to a large contribution 
to the temperature gradient by changing from a-0.5 to a-1.25 
if that happens in an over-adiabatic layer. Moreover, we must 
take into consideration the fact that the Kurucz grid is sparser 
than the NEMO one, and the derivatives with respect to T e s 
and log g involved in </> x and A x are hence determined with a 
larger uncertainty. These facts can be seen in Figs. [T7l and [TBI 
which are similar to Figs. [15] [16] but for a slightly different 
1.8 M evolutionary state. These stellar models are only 50 K 
cooler than the previous ones, and their gravity is different by 
A log g = 0.014. Now, the behavior of A x /A y and X - y for the 
MLT(a=0.5)+Kurucz-MLT(a=1.25) model has changed, and 
the differences with FST models have significantly increased. 

As the amplitude ratios are concerned, both the atmosphere 
models and non-adiabatic results can play a significant role, 
mainly for higher degree modes. Information on the flux and 



Fig. 15. Theoretical monochromatic phase differences <p u ,v,b,y ~ 
(p y as a function of wavelength for 1.8 M Q models with 
r e ff=6800 K and logg=3.816. Each panel corresponds to a 
different degree i for frequencies close to the fundamen- 
tal mode. Solid line: FST interior, FST atmosphere struc- 
ture and color transformation. Dashed line: MLT(ar=0.5) in- 
terior, MLT(ar=0.5) atmosphere structure and color transfor- 
mation, dotted line: MLT(q'=0.5) interior, MLT(q'=0.5) atmo- 
sphere structure and color transformation and limb-darkening 
coefficients derived from FST atmosphere models (labeled by: 
FST LDC). Dot-dashed line: MLT(a:=0.5) interior, Kurucz 
MLT(a=1.25) atmosphere structure and color transformation. 

limb-darkening could thus be extracted from the amplitude ra- 
tios if the non-adiabatic predictions are known with a suffi- 
ciently high degree of confidence. 

4. Application to the 5 Scuti star AN Lyn 

AN Lyn is a medium amplitude 5 Scuti star discovered by 
Yamasaki et alj (Il98ll) . Subse quent investigations revea led that 
it is a multi-periodic 5 Scuti (Rodriguez et al. 1997b) with a 



peculiar light curve (Rod riguez et al.lll997al) . iRodriguez et al 
dl997bl) identified three independent frequencies: v\ = 
10.1756 c/d, v 2 = 18.1309 c/d and v 3 = 9. 5598 c/d, an d 



two of them (vi and V2) have been confirmed bv lZhoul(l2002h . 
Moreover, the amplitude of the dominant frequen cy (vi ) seem s 
to show variations with a time-scale of years dZhoull2002l) . 
Precise amplitudes and phases in Stromgre n uvby photome- 
try were derived by Rod riguez et al. ( 1997b), we will consider 
these observational results to apply our theoretical models. 

The Stromgren i ndices, together with the c alibration soft- 
ware TempLoGv2 (IStiitz & Nendwichl |2002J), provide the 
following global parameters: r e ff=7150 K, log,g=3.65 and 
[M/H]=-0.17, and taking into account six different calibra- 
tions: r eff =7100 ± 150 K, logg = 3.5 + 0.35 and [M/H]=- 
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Fig. 16. As Fig. [15] but for monochromatic amplitude ratio Fig. 18. As Fig. [16] but for models 1.8 M Q models with 



r eff =6750 K and logg=3.802. 
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Fig. 17. As Fig. [15] but for models 1.8 M models with 
r eff =6750 K and logg=3.802. 



0.17+0.03. Even if the metallicity resulting from this calibra- 
tion is sub-solar, we shall do the computations with classical 
solar composition, and NEMO atmospheres for [M/H]=0 and 
microturbulence velocity £ t = 2 kms" 1 . There are several rea- 
sons for these choices: first of all, the value of metallicity pro- 
vided by the photom etric calibration is quite uncertain. Second, 
Barban et alj (120031) computed the limb-darkening coefficients 
only for that metallicity and £ t ; and finally, the goal of this sec- 



tion is not a precise fitting of the AN Lyn observations but to 
analyse how the choice of different treatments of the convective 
outer layers may affect the mode identification. 




350 400 450 500 550 350 400 450 500 550 
Wavelength (nm) Wavelength (nm) 

Fig. 19. Photometric amplitude ratios A u ^^^jA y (lower panels), 
and photometric phase differences <p u ,v,b,y - <Py (upper panels) 
for different models of AN Lyn. Lines refers to theoretical pre- 
dictions for I — (solid line), I — 1 (dashed line), 1-2 
(dot-dashed line) and I — 3 (dotted line), and the error bars 
represent the observations for v\ . 
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We have computed a small grid of stellar models for Z = 
0.02, X = 0.7 chemical composition, and masses from 1.70 to 
2.20 M . The best fit of r eff (7166 K), logg (3.76), main fre- 
quency and photometric features (amplitude ratios and phases) 
is obtained for M=2.0 M at t = 1.10 9 yr, X c = 0.145, a mass 
of the convective core M cc /M = 0.18, and a density contrast 
p c /p ~ 860. As in previous papers concerning AN Lyn, we 
identify vi as a n — 2 radial mode. Its second frequency how- 
ever is affected by too large errors to be surely identified. For 
this 2 M Q model, this frequency could correspond to a £ = 3 
n = 2 mode. 

The four panels of Fig. [19] show for the different (uvby) 
bands the amplitude ratios and phases, which were computed 
for complete MLT and FST models, and transformed to the ob- 
servational plane by using the limb-darkening coefficients and 
integrals for the corresponding model atmosphere. It have also 
checked the effect of combining a kind of atmosphere (MLT) 
in the stellar model with the color and limb-darkening transfor- 
mations derived from FST model atmosphere. As it was already 
pointed out, the main effect is due to the temperature gradient 
in the external layers, and the effect of using FST or MLT based 
limb-darkening coefficients is negligible. 

We have also computed MLT(a=1.25) models, using con- 
sistently the Kurucz model atmosphere grids and the corre- 
sponding limb-darkening coefficients. The temperature gradi- 
ent derived from this models are unable to predict the expected 
amplitude ratios, especially in the u band. 

Concerning FST or MLT (a=0.5), the differences are very 
small and not observationally disentangled. This result is not 
surprising, given the small differences in the theoretical non- 
adiabatic quantities shown in Figs.|9]and [l0]for a 2 M star. 



5. Conclusions 

Mode identification methods for 6 Scuti pulsations are based 
on the non-adiabatic quantities fj and tfrj which are highly sen- 
sitive to the characteristics of the surface convective zone. The 
observed phase-lag between light variation and the displace- 
ment originates in two different regions: the Hell partial ion- 
ization zone, where the k mechanism drives oscillations, and 
the HI and Hel partial ionization region, which coincides with 
a convective zone (HCZ). While the contribution from the Hell 
region mainly depends on its depth, and therefore on the evo- 
lutionary state, that of HCZ is very sensitive to the temperature 
gradient as well and, therefore, to the convection treatment. In 
this paper we have studied the role of convection in the ex- 
ternal layers by comparing low efficiency convection models: 
FST and MLT(ar=0.5). The novelty of our analysis is that, for 
the first time, the model atmosphere used in the pulsation anal- 
ysis were computed with the same treatment of convection as 
in the internal structure. 

This study has shown that: 

- Even if FST and MLT(a;=0.5) model atmosph ere provide 



ers of FST and MLT(ff=0.5) over-adiabatic regions, show- 
ing differences of ~ 20° in the phase lag. 

- Differences increase as decreases. This is due to the 
different convective fluxes in FST and MLT treatments. 

- For low efficiency of convection, the interaction pulsation- 
convection is negligible. So, if convection efficiency is low 
enough, we find for the "frozen convection" approximation 
phase-lag values close to results obtained taking into ac- 
count the convection-pulsation interaction. Given that FST 
is less efficient than MLT, i/^r for FST remains higher than 
100° for lower r eff than MLT(ff=0.5). 

- Differences between FST and MLT(ar=0.5) also decrease 
as stellar mass increases. With increasing stellar mass, the 
HCZ becomes more and more shallow and, therefore, con- 
vection less and less efficient. 

- The FST results cannot be reproduced with a single a 
value, /t and iffj from FST are bracketed by MLT(a=0.1) 
and MLT(ar=0.3), but depending on the frequencies the 
FST results are closer to MLT(ar=0.1) or to a=0.3 
ones. Furthermore, for higher r e jf the FST behavior is 
closer to that of MLT(ff=0.1), but becomes closer to the 
MLT(a;=0.3) one as T e ff decreases. 

- The use of Kurucz model atmospheres with ar=1.25 in stel- 
lar models with <?=0.5 have shown that the differences due 
to the very external layers can introduce an uncertainty 
in the phase lag of the order of 40° when compared with 
"complete" MLT(a=0.5) models. 

- As for iff j, photometric phase differences in the Stromgren 
system are mainly sensitive to the temperature gradient in 
the HCZ, and are only slighly affected by the color trans- 
formation tables and limb-darkening functions provided by 
different model atmosphere. The amplitude ratios, on the 
contrary, are affected by both the thermal structure and the 
color transformations. The u and v bands are in general the 
most affected: differences in <p u - <p v of the order of 3-4° 
may be expected for {=0 or 2; and differences of 10% to 
20% for the amplitude ratio values. 

- Differences between FST and MLT(ar=0.5) monochromatic 
phase differences and amplitude ratios are also frequency 
and t dependent. For frequencies close to the fundamen- 
tal mode, the sensitivity to the convection treatment does 
not prevent the identification of the mode degree Z from the 
classical amplitude ratio vs. phase differences diagram. But 
as frequency increases, ( — 0, 1 or 2 mode identification be- 
comes difficult due to the uncertainty on FST/MLT convec- 
tion treatment. Adding radial velocity data could help, nev- 
ertheless, on discriminating between FST and MLT(ar=0.5) 
approaches. 

- The application of our study to the 8 Scuti star AN Lyn pro- 
vided results in good agr eement with previous publications 



(Rod riguez et al 



similar stellar spectra and limb-darkening (IHeiter et al 



20021: iBarban et alj2003b . the non-adiabatic observables f\ 



and clearly reflect the thermal structure of the inner lay- 



1997a). Moreover, given the high mass 
derived for this star (~ 2.0 M ) the effect of convection 
treatment is small, and it is not possible to discriminate be- 
tween FST and MLT(ar=0.5) on the base of the available 
observations. 



In summary, the theoretical photometric amplitude ratios 
and phase differences between different photometric passbands 
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are directly related to the non-adiabatic results (see Eq. [TJ. 
Moreover, different modes probe different layers of the star. 
Hence, the temperature gradient in the superficial convection 
zone of intermediate mass stars can be constrained by compar- 
ing these amplitude ratios and phase differences with observa- 
tions. Such information cannot be obtained from the spectrum 
alone because spectral lines probe only region of small optical 
depth. 
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